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o ■ 

^vQ , We investigate two- and three-particle intensity correlation functions of pious in rela- 

tivistic heavy ion collisions for different colliding energies. Based on three models of particle 
production, we analyze the degree to which the pion sources are chaotic in the SPS S-|-Pb, 
Pb-|-Pb and the RHIC Au-I-Au collisions. The "chaoticity" , A, of the two-particle correla- 
tion functions is corrected for long-lived resonance decays. The effect of the partial Coulomb 
correction is also examined. Although the partially coherent model gives a result which is 
consistent with that of STAR, the chaotic fraction does not exhibit clear multiplicity de- 
pendence if we take into account both the corrected chaoticity and the weight factor of the 
three-pion correlation function. The result of the partially-multicoherent model indicates an 
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• increasing number of coherent sources in higher multiplicity events. 

. §1. Introduction 

o ■ 

' Pion interferometry has been regarded as an indispensable tool in relativistic 

heavy ion physics. Two-particle intensity interferometry can be used to determine 
the sizes of the collision system. This fact is known as the Hanbury Brown-Twiss 
(HBT) effect. For this reason, it has been used for exploring the space-time evolution 
^ . of hot, dense matter created in heavy ion collisions.^-' In particular, the most recent 

pH I experiment at the Relativistic Heavy Ion Collider (RHIC) obtained an interesting 

result which is referred to the "HBT puzzle" .^^'^^ Hydrodynamical models have 
failed to reproduce the experimental results of two-particle correlation functions. 
r> I Hydrodynamical models are based on the assumption of local thermal equilibrium. 

■ For the early stage of the space-time evolution, the validity of this assumption is 

indirectly verified by the observation of the large elliptic flow at RHIC.^^ However, 
equilibration in the final hadronic stage is still ambiguous. Although an exponential 
particle spectrum is expected for a thermal source, it does not require a system that 
has reached local thermal equlibrium.^) 

The two-particle intensity correlation function C2{pi,P2) of identical particles 
provides information not only on source sizes but also on the state of the source 
through the chaoticity, A = C2(p,p) — 1. The HBT effect does not exist if the source 
is coherent. The chaoticity A is unity for a perfectly chaotic source and for a 
coherent source, which does not induce the HBT effect. In multi-particle production 
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phenomena of relativistic heavy ion colhsions, the thermahzed source can be chaotic 
and non-thermal components can be coherent. For example, if disoriented chiral 
condensate domains are created, such domains can decay into coherent pions. Thus 
the final state pions may carry information regarding how partons hadronize if the 
deconfined phase has been created. Therefore, the chaoticity A is a very important 
quantity to investigate the final state in heavy ion collisions. 

However, A cannot be regarded as the true chaoticity in real (experimental) 
situations, because many other effects, such as long-lived resonance decay contri- 
butions, ^^"^^ Coulomb repulsions, and particle contaminations^''-* affect the chaotic- 
ity A. As an alternative tool, the three-particle correlation function has been pro- 
posed. ^^^"^^^ Three-particle correlations are more useful for this purpose, because 
long-lived resonances do not affect the normalized three-pion correlator, 

^3(Pl,P2,P3) 

_ [C3(Pl,P2,P3) - 1] - [g2(Pl,P2) - 1] - [C2{P2,P3) - 1] - [C2{P3,Pl) ' 1] 
V[C72(Pi,P2) - mC2{P2,P3) - l][C2(P3,Pl) " 1] 

(1-1) 

with C3{pi,p2,Ps) being the three-particle correlation function. The correlation 
functions are defined as 

and 

Wl{p-i_)Wl{p2)Wl{p^) 

with Wn{pi,-" iPn) being the n particle distribution. The index of the source 
chaoticity in the three-pion correlator is the weight factor 

uj = r3{p,p,p)/2 (14) 

which is unity for a chaotic source. Due to insufficient statistics, the three-pion 
correlator measured in experiments to this time is slightly different from Eq. 

^(Q)= [g3(g3) - 1] - [C2(Ql2) - 1] - [C2(Q23) - 1] - [C2(Q3l) - 1] 
V[C2{Ql2) - 1][C2{Q23) - l][C2{Q3l) " 1] 



where Qij = ^ —(j>i — pj)'^ and Qs = \/ Q12 + Q23 + Qsi- The weight factor is 
defined as w = r3(0)/2, so that it is expected to be the same as the definition 
Eq. (II31). 

As shown in a previous work,^^) more precise information concerning the source 
can be extracted using model analyses combining two- and three-particle correlations. 
In this paper, we analyze the two-pion and three-pion correlations and investigate 
the chaoticity of the pion sources. Firstly, assuming that the main background con- 
tribution is a long-lived resonance decay and other effects are successfully removed 
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in the experimental data, we make a correction to account for the long-hved reso- 
nance decays to the chaoticity of the two-pion correlation functions, with the help 
of a statistical model. Secondly, we extract the weight factor to through the simul- 
taneous construction of C2 and C3. Finally, we carry out model analyses using the 
"true" chaoticity A*'^'^'^ after applying the resonance correction and the weight factor 
uj. In addition to the analysis of Au-|-Au collisions at the RHIC given in a previous 
paper, ^''^ we give analyses on lower colliding energy collisions at the CERN Super 
Proton Synchrotron (SPS) based on the same method. We treat data from 200A 
GeV (laboratory system) S-|-Pb collisions measured by the NA44 coUaboration^^^"^*') 
and data from 158^ GeV (lab. sys) Pb-|-Pb collisions measured by the NA44 collab- 
oration^^)' and the WA98 collaboration.^^)' On the basis of the results obtained 
from these experimental data, we investigate the multiplicity dependence of the ex- 
tracted model parameters. 

In the next section, we explain the three models used in this paper. The cor- 
rection to account for the long-lived resonance decays is discussed §3. In §4, we will 
give a combined analysis of the 27r and Svr correlation functions, with the goal of 
extracting the weight factor lo. The results of the model analyses are given in §5. 
The paper is summarized in §6. 

§2. Model description 

We use three kinds of models in this analysis. In all the models, A*"^"*^ and uj are 
used as inputs to fix the model parameters. The output quantities are the chaotic 
fraction and the mean number of coherent sources. Here we do not discuss the origin 
of the coherences in these models. Figure ^ presents a schematic depiction of the 
models. 

One is a partially coherent model^^^ whose only parameter is the chaotic fraction, 
ei, defined as the ratio of the number of particles coming from the chaotic source to 
the total number of particles. In this model, pions are emitted from a mixture of a 
chaotic and a coherent source. In general, we need to fix the source function in a 
manner that depends on the momentum. In the present analysis, fortunately, the 
true chaoticity A'''^'^'^ and the weight factor oj are given by the correlation functions at 
vanishing relative momenta, and hence they can be expressed in terms of the chaotic 




Fig. 1. Schematic depiction of Models I-III. The shaded area denote the chaotic sources. The 
unshaded area and the small circles represent the coherent sources, (a) Model I, (b) Model II, 
(c) Model III. 
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fraction. The relations among £i, A and uj are 

A^™^ = ei(2-£i), - = V^(^^- (2-1) 

We refer to this model as Model I. 

The second model, referred to as Model II, is a multicoherent source model. 
The parameter in this model is the mean number of coherent sources, an, which 
obeys the Poisson distribution. Because small coherent sources are independent, a 
chaotic source is realized as a cluster of an infinite number of coherent sources. The 
parameter an is thus related to A*"^"^ and u as 

x^true _ "II ^ _ 1 '^<A + 2«ii + 3 fan + 1 ^2.2) 
an + 1 ' 2 + Sail + 1 V an ' 

As shown in Fig. 2 of Ref. 14), w as a function of an has a minimum value o^min — 
0.82. Hence, there is no corresponding an for u smaller than tJmin- 

As each of these models possesses only a single parameter, the parameters cal- 
culated from A*'^"^ and u should give the same value if all corrections are correctly 
made. 

Model III is a "partially multicoherent" source model.^^^ This model is a mixture 

of Model I and Model II, and its parameters are the chaotic fraction eni and the 
mean number of coherent sources am, which are related to A*''"'^ and u) as 

A*™^ = ^, (2-3) 

am + (1 - eiii)^ 



^ _ 2ai^n + 2am(l - em? + 3(1 - £m)-^(l - 2£m) / am + (1 - giii)^ 
2[a2jj + 3am(l-eiii)2 + (l-em)3] V «i" ' 

In this model, there are two parameters that correspond to two experimentally mea- 
surable quantities, A*"""^ and uj. Below, we investigate the allowed parameter regions 
for given sets of A*''"^ and u. 

§3. Extracting A*'^"'' from two-pion correlations 



In relativistic heavy ion collisions, a non-negligible fraction of pions comes from 
the decay of long-lived resonances. Recent analyses based on statistical models show 
that hadrons are chemically frozen near the critical temperature Tg.^^-* Short-lived 
resonances, such as p and Z\, decay before hadrons reach the kinetic freeze-out or 
soon after the freeze-out, but long-lived resonances, such as hyperons, can decay long 
after the kinetic freeze-out of pions. In the two-pion correlation function, these long- 
lived resonances appear as a spike near (? ~ 0, whose width is too small to be resolved 
with the current experimental resolutions, Aq ~ 5 — 10 MeV. Thus, chaoticities A^'^p 
measured experimentally are smaller than the true chaoticities due to long-lived 
resonance decays.^)' ^^'^^^ Following Ref. 26), we take into account resonances up to 
Z'*(1385). We treat resonances whose widths are less than 5 MeV as long-lived ones, 
i.e., -fC^, ri,ri',(j),A,E and ^ are considered long-lived resonances in the calculation. 
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Table I. Thermodynamic parameters obtained from the particle ratio in nucleus-nucleus collisions 
at SPS. Here, h~ represents the negatively charged hadrons. 



System 


Ratio 


Data [Reference] 


T [MeV] 


MB [MeV] 




SPS S+Pb 


P/P 
P/tv~ 
{K- + K+)/2K° 
K+/K- 

kVa 

A/ A 

S+/A 


0.18±0.03"»)'"'') 
0.12±0.02^°) 
0.024±0.009^*''^''' 
1.07±0.03^^^ 
1.67±0.153^> 
1.4±0.1^^' 
6.4±0.4^^' 

0.20±0.01^^' 

0.21±0.02='2) 


173 ± 2 


196 ± 2 


36/7 


SPS Pb+Pb 


p/p 


0.085±0.009^^> 


161 ± 4 


223 ± 7 


44/9 






0.125±0.019^''^ 








K°/h- 


0.123±0.020^^^ 










A/h- 
A/K° 


0.077±0.011^^^ 










0.63±0.08^^^ 










A/A 


0.131±0.017''*"''^ 










E-/A 


0.110±0.010=^5) 










S/A 


0.188±0.039^^' 










{S + S)/{A + A) 


0.13±0.03^''^ 












0.232±0.033^®^ 










K+/K- 


1.85±0.093*' 









Though u mesons have an intermediate width, it is known that they distort the 
shape of the correlation function at low q but do not reduce the chaoticity.'''^ 

For a chaotic source, the reduced chaoticity A*^^ is given in terms of the ratio of 
the number of pions from the long-lived resonances to the total number of pions as 

VAa = i-f, (3.1) 

where N.,^ is the total number of emitted pions, and is the number of pions from 
the decay of long-lived resonances.^) We calculate this ratio with the help of the 
statistical model. For midrapidity, the particle ratio can be written as the ratio of 
the number densities, i.e., 

Ni n? 



N, n. 



o> (3-2) 



where 

^ 27r2 



poo 

I dpp'^f{E,T,iJ,B,IJ'S,IJ'i3), (3-3) 
^0 



with f{E,T,fXj,---) being the equilibrium distribution function and Qi being the 
number of the degree of freedom of particle z.^^^ For simplicity, wc fix the chemical 
potential of the third component of the isospin as /ij^ = 0. Then, the thermodynamic 
parameters to be determined are the temperature T and the baryon number chemical 
potential fiB, because the strangeness chemical potential fis is determined from the 
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Table II. Thermodynamic parameters obtained from the particle ratio in Au+Au collisions at 
RHIC. 



System 


Ratio 


Data 


r [MeV] 


liB [MeV] x"/Nrio{ 


RHIC Au+Au 


P/P 


0.71±0.05^^' 


158 ± 9 


36 ± 6 2.4/5 




p/tt" 


0.072±0.014*"''^^' 








K~/tv- 


0.146±0.02*2) 








K*°/h~ 


0.042±0.014''3' 










0.92±0.14''^) 








A/A 


0.71±0.05''^' 










0.83±0.09''^^ 







strangeness neutrality condition. Results for several collision systems obtained from 
the fit are shown in Tables HI and ITll*^ 

From Tables and ^ we find that the quality of the fit becomes better as 
the collsion energy increases. It is reasonable to assume that the system has reached 
thermal equilibrium and that the pion source becomes chaotic as the collision energy 
increases. Below we see whether this naive assumption is valid. 

Since this particle ratio is obtained from the particle numbers integrated with 
the particle momenta, A'^''^ is calculated from Eqs. H3-l|) - H3-3|) using the integrated 
particle numbers. The measured chaoticity A'^^p, however, depends on the average 
momentum of pion pairs. For example, STAR has measured X^^^ for three bins of 
the transverse momentum (/ct).^^^ In this paper, we assume, for simplicity, that 
the true chaoticity A*"^"*^ does not depend on the particle momentum. Therefore, 
we need to average X^^^ over the momenta in order to evaluate A*'^"^. Assuming 
the kt dependence of the measured A^^p is dominated by long-lived resonances [i.e., 
Eq. (|3-lj) ] , we obtain the averaged chaoticity A™'' as 



n „ 



ktdkt 



dN 



2 



—pxn ,-1 ^i-thbin \ktdkt _ 

^ = ; (3-4) 

ktdkt 



E 

1=1 



i-th bin 



ktdkt 



with A^^'' being the measured chaoticity in the i-th kt bin. [See the appendix of 
Ref. 17) for the derivation of Eq. ((231) •] The transverse momentum distribution 
dN/ ktdkt, in the above equation is taken from experimental results. For S+Pb col- 
lisions, in which three-particle correlation data are available only for minimum-bias 
data, we calculate '}^^^ by simply averaging the A'^^p with different multiplicities. 
This should not affect our conclusion, because the data exhibit little multiplicity 
dependence. For A'^^p, it is known that the value of 6*2(9) — 1 at q = depends 
on the dimension of the projection onto q-space. In general, the value of A'^^p ob- 
tained from the 1-dimensional Gaussian fitting differs from that obtained from the 
3-dimensional Gaussian fitting because of the projection, experimental resolution, 

*-* Though experimental data can contain contributions from coherent sources, we assume that 
particle ratios are not affected by the existence of coherent sources. 
**' Note that this is also assumed in Model I-III. 
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Table III. Summary of A"'"'^ 



System 


Experiment [Reference] 




^truc 


A true 


SPS S+Pb 


NA44, min.bias^^' 


0.585 ± 0.06 


0.94 ± 0.06 


O.TA"^"" 


SPS Pb+Pb 


NA44 Central^i) 


0.55 ± 0.03 


0.98 ± 0.03 


0.8A"^"'= 


SPS Pb+Pb 


WA98 CentraP"*) 


0.58 ± 0.04 


1.03 ± 0.04 


0.8A"^"'= 


RHIC Au+Au 


STAR Central*''' 


0.57 ± 0.06 


0.93 ± 0.08 


0.8A''"'= 



and other effects, though these should be the same for an ideal measurement. In this 
paper, we use the value extracted from the 3-dimensional Gaussian fitting, 



C2*(q) - 1 + ^''^^ exp (--Rside^sido - -^out^out " -^longO'long) • (3-5) 



as A^^P. 

The true chaoticity is then given by A*"^"*^ = 'X'^^ /X^^. Results for A*'^'^'^ obtained 
from various systems are summarized in Table ITTTI The error on A*"^"*^ is the sum of 
the experimental one on A''^'', calculated from the errors on A^^'' and dN/kfdkt, and 
the errors propagated from the fit of the thermodynamic quatities at the 1-a level, 
shown in Tables HI and ITTl 

From Table |llTl we see that the value of A are not so different and become 
close to unity in all systems. However, it should be noted that there may be an over- 
estimation in A^^P due to an overcorrection for the Coulomb interaction. Because the 
two-pion correlation function is affected by the Coulomb interaction between the two 
detected pions, there are many issues concerning how the Coulomb interaction can 
be subtracted from the observed correlation function. ^^^^^'^^ Recently, it has been 
shown that a new procedure called the partial Coulomb correction, leads to signif- 
icant correction to the source sizes. Such corrections are known to partially 
resolve the "HBT puzzle", that -Rout becomes smaller than iigide; though this correc- 
tion does not completely resolve this puzzle. In this paper, we stress that the partial 
Coulomb correction affects not only the source size but also the observed chaoticity, 
^exp_ For example, CERES reports that the kt dependent correction to A'^^p reaches 
15-40 This is a significant correction. Unfortunately, partially corrected data 

for the collisions that we treat in this paper (S+Pb, Pb+Pb at the SPS and Au+Au 
at 130^ GeV) are not available. Ideally, these data should be treated within a 
theoretical approach with a Coulomb correction, but this is beyond the scope of 
this paper. In the further analyses given in ^ we also present the results obtained 
using "partially corrected" A'^^p in order to see how the results change when A°^p is 
reduced. Following the report from CERES, ^'^^ according to which the correction is 
larger at smaller kt , and considering that the two-particle data used here are those of 
the lowest momentum bin, we simply multiplied a correction factor of 0.8 except in 
the S+Pb case, where the correction factor is set to 0.7, because only S+Pb data are 
corrected by the Gamow factor. We denote this "partially Coulomb corrected" A*'''^'' 
by A^r • 
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§4. Extraction of the weight factor u: 

4.1. Procedure 

In the previous section, we extracted A*''™ from the two-particle correlation data. 
Next, we must determine the weight factor uj in order to analyze the degree to which 
the pion sources are chaotic by using the models, especially Model III [Eq. (|2-3() ]. 
To obtain uj = r3(0)/2, we must extrapolate C2 and C3 to q = 0. The method of 
extrapolation is described in a previous paper^'') in detail. Here we briefly review 
the procedure. 

In experimental papers, uj is extracted by simply averaging r3{Q-^) (NA44,^'^)'^^^ 
WAOS^^-*), in which there is little dependence on Q3, or using quadratic and quartic 
fits (STAR^^^). A shortcoming of these methods is that the Q3 dependence of 
may be more complex than quadratic or quartic. Such Q3 dependences result 
from both an asymmetry of the source and coherence.*^ In this paper, we assume 
observed Q3 dependences are due to coherence. This is an assumption, but it is 
plausible, because the asymmetry of the source causes a Q3 dependence of that is 
somewhat different from that in the observed r^ data.^^^ We reproduce C2(Qinv) and 
CsiQs) using a common source function with a set of parameters which is determined 
by minimizing with respect to the experimental data. Then, we evaluate 6*2(0) 
and C3(0). The quantity r3(Q3) is also calculated for a consistency check. For 
simplicity, we use a spherically symmetric Fourier-transformed source function with 
simulateneous emission, Fij = fij{\qij\)e^^^'-~^i^^° , in which the exponential term 
corresponds to emission at a constant time to. The assumption of simulateneous 
emission should be a good approximation, because the experimental data suggest 
emissions of short duration throughout the broad range of colliding energies. Since 
the finite emission time duration is not related to A*™*^ but, rather, to the width of 
the outward correlation functions, it does not affect our results below. For the spatial 
part, fijilQijl), we try the following three kinds of source functions: 



h,ij{\q,^\)=exp{-R^\q,//2), (4-1) 
/2,.,(|q,,|)=exp(-fi|q,^.|/2), (4-2) 

/3,,(|q,,|) = ^ ^ (4-3) 
cosh{R\qij\) 

Here, q^j = Pi — Pj and i? is a size parameter which is to be determined by the 
fitting. The third function, Eq. (|4-3|) . is chosen so as to be quadratic at small |q| 
and exponential at large |g|. 

The two- and three-particle correlation functions are then calculated using the 
relations 



*' In Ref. 52), it is pointed out that there is a possibility that the standard projection method 
adopted here may lead to an artificial momentum dependence of the projected correlators. 
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Table IV. Results of the fitting to C2 and C3 



system 




R [fm] 


'^inv 




xVdof 




SPS S+Pb 


fi 


4.85±0.31 


0.49±0.04 


0.34±0.04 


2.3/7 


0.33±0.38 


(NA44)i«)'20) 


/2 


7.55±0.84 


0.79±0.10 


0.61±0.09 


3.9/7 


0.48±0.45 




/3 


8.99±0.80 


0.55±0.06 


0.40±0.05 


1.9/7 


0.40±0.44 


SPS Pb+Pb 


/l 


7.37±0.61 


0.54±0.05 


0.49±0.07 


1.2/6 


1.06±0.59 


(NA44)2i)'22) 


/2 


lf.4±1.5 


0.83±0.11 


0.85±0.16 


2.7/6 


1.16±0.69 




/3 


i3.8±1.5 


0.6f±0.07 


0.58±0.10 


0.9/6 


1.15±0.67 


SPS Pb+Pb 


/l 


7.73±0.10 


0.36±0.01 


0.29±0.0f 


113/17 


0.81±0.12 


(WA98)24) 


/2 


14.2±0.29 


0.75±0.02 


0.62±0.02 


14/17 


0.68±0.12 




/3 


15.8±0.28 


0.47±0.01 


0.39±0.0f 


17/17 


0.78±0.14 


RHIC Au+Au 


/l 


7.0±0.07 


0.54±0.01 


0.48±0.01 


110/30 


0.958±0.09 


(STAR)4'')'5i) 


/2 


i4.4±0.2 


1.18±0.03 


1.08±0.03 


79.7/30 


0.736±0.09 




/3 


i5.2±0.2 


0.71±0.01 


0.64±0.02 


15.8/30 


0.872±0.097 



/II/22 

C3(P1,P2,P3) = 1 + - f E + 2-3^^^) . (4-5) 

\ ^ iiilji JUJ22J33 I 

where Ainv, and 1^3 are phenomenological adjustable parameters accounting 
for the non-trivial coherence effect. The summation ^(jj) here runs over {i,j) = 
(1, 2), (2, 3), (3, 1). The quantities Amv, and are unity in the case of a fully 
chaotic source. We can set 1^3 = 1 for a description of Cs{Qs) at small Qs-^^^ Then, 
Ainv and v are determined by the fitting, like R. We stress that the fit is 
carried out simultaneously for the two- and three-particle correlation data. 

The results of the fittings to the experimental data and the resultant uj are 
listed in Table HB 

4.2. SPS S+Pb 

Comparisons of the fitting results and the experimental data are displayed in 
Figs. I2ta), Ota) and 01 a). The curves in Figs. HI a) andOl^a) are fitted to the data 
using three kind of source functions: Gaussian [Eq. ()4-l() ]. exponential [Eq. ()4-2() ] 
and cosh [Eq. ()4-3() ]. The fitting range is < Qij^Qs < 60 MeV, which is adjusted 
to the available Q3 range for r3((53). In the case of C2, we see that the experimental 
data are fit well by a quardratic function at small Qmv and an exponential function 
at large Qmv For this reason, the values in the exponential case and the Gaussian 
case are larger than those in the cosh case. A similar tendency is also seen in the case 
of C3, though the exponential case is not excluded by the data point at the smallest 
Qinv The value of r^ obtained from the fitted C2 and C3 do not differ greatly. With 
a simple naked-eye extrapolation in Fig. Ufa), it appears that uj = r3(0)/2 becomes 
~0.4 for all source functions. Taking the smallest x^ value, we adopt the cosh case 
and obtain uj = 0.40ib0.44. Note that, by the definition of r3((53), where C2{Qij) — 1 
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rnt^h v^/rlnf-1 Q/7 

X=0.55 

^\ exp, Y^/dof=3.9/7 

\ ?.=0.79 

~~r<^ gauss, Y7dof=2.3/7 
NA44 3i"^ Min.bias 

1 1 1 1 1 1 1 = 


(h'\ cosh y^/dof=0 9/6 

, X=0.6^ 

\ exp, x^/dof=2.7/6 

■O, ?.=0.83 

t\ gauss, Y^/dof=1 .2/6 
%i ?L=0.53 

\^ NA44jt" ^ 


(c) cosh, Y^/clof=1 6.8/1 7 

. X=0.47 

\ exp, x^/clof=14.2/17 --- 

\ gauss, x^/clof=113/17 
~nA^^ X=0.36 

— WA98jr" ^ 


■J — 1 1 1 1 1 1 1 1 

\ (A\ cosh, x^/dof=1 5.8/30 

■ \^°^ X=0.7^ 

X\ exp, x^/dof=79.7/30 --- 

■ \ X=1.18 

■■^ gauss, x^/clof= 11 0/30 
V ?L=0.54 

^ STAR 7t" • 



10 20 30 40 50 60 70 80 
Qinv [MeV] 



10 20 30 40 50 60 70 80 90 
Qinv [MeV] 



Fig. 2. Two-pion correlation function C2(Qinv) in the (a) S+Pb collisions at the SPS, (b) Pb+Pb 
collisions at the SPS, measured by NA44, (c) Pb+Pb collisions at the SPS, measured by WA98, 
and (d) Au+Au collisions at the RHIC. The lines represent our results for the fits of each source 
function (see text). Filled circles represent the experimental results (SPS S+Pb by NA44,'^^-' 
SPS Pb+Pb by NA442i)and WAgS,^*' and RHIC Au+Au by STAR''"''). 



is in the denominator, there exists a large uncertainty in the calculated ?'3(<33) value, 
because of the errors on Ainv In Fig-S^a), we plot the error band for the cosh case. 
The uncertainty 0.44 associated with lo is also obtained by extrapolation of the error 
band to Q3 = 0. 

4.3. SPS Pb+Pb 

For Pb+Pb collisions at 158A GeV, we compared our results to data measured 
by both NA44 and WA98. Figures Efb), Ofb) andlljb) plot the results of the fit 
to the NA44 data and Figs. I^c), Ofc) andH^c) plot those for the WA98 data. The 
fitting ranges are adjusted to the available r3{Q3) data range, as in the S+Pb case, 
and we use the data for Qij, Q3 < 60 MeV. 

Comparing our fit to the NA44 data, we find that all source functions seem 
to give nice descriptions of the data [see Figs. H^b) andOl^b)]. The statistics are 
still insufficient, especially at low Q^, in the three-particle correlation function to 
discriminate the best source function. The weight factor uj is larger than unity at 
the best fit value. This is associated with the large errorbars and consistent with 
the WA98 case within the errorbars. For further analysis (given below), we adopt 
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Fig. 3. Three-pion correlation function C3{Q3) in the various 
symbols is the same as in Fig. |21 The experimental results 
Refs. 20), 22), 24), and 51), respectively. 



collisions. The identification of the 
presented in (a)-(d) are taken from 



the cosh case because it gives the best value. On the other hand, the WA98 data 
exhibit better statistics [see Figs. [5Ic) and|n[c)]. This excludes the Gaussian case 
in the fit to C2 and C3. For the values, the exponential case has the best value, 
X^/doi = 14/17. This tendency seems to be different from that in the NA44 case. 
In Ref. 24), the three-particle correlation function is fit well by a double-exponential 
correlation function. Although the exponential case seems to be the best of the three, 
faiQs) from the exponential source function deviates from the experimental result 
[Fig. Etc)]. This should not be regarded as a serious problem, however, because 
errors on Ajnv and v lead to an uncertainty on r^^Q^), as shown in the SPS S-|-Pb 
case [Fig. Efc)]. A likely reason for this deviation is that C3{Qs) is smaller than 
that obtained from the experimental data in the exponential case at low Q^. Naive 
extrapolation by naked eye in Fig.^c) again suggests that the cosh or Gaussian case 
is better, but this may be a coincidence in the Gaussian case, because both C2 and 
C3 deviate from the experimental data at low relative momenta. Hence, we adopt 
the cosh the result for w. Despite the different behavior of C3{Q3), the result 

for uj obtained from the WA98 data is consistent with the that obtained from the 
NA44 data. 
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Fig. 4. The normalized three-pion correlator r3{Qz)/2. The identification of the symbols is the 
same as in Fig. |5| The experimental results presented in (a)-(d) are taken from Refs. 20), 22), 
24), and 51), respectively. The dot-dashed line in (a) represents the uncertainty propagated 
from the fitting parameters for the "cosh" case (see text). 



4.4. RHIC Au+Au 

Results presented in this subsection are the same as those presented in Ref. 17). 
Figures [2td), Ofd) and 0fd) display the correlation functions and the three-pion 
correlator for the Au+Au collisions at the RHIC energy. The fitting range is < 
Qinv < 90 MeV for C2 and < Q3 < 100 MeV for C3. The high statistics of the 
data allows us to discriminate the source function. In the cosh case, the value of 
is much smaller than other two cases. Finally, we obtain u: = 0.872 it 0.097. 

§5. Results and discussions 

In ^ we extracted A*''"'' from the experimental data, A'^^p, with the help of a 
statistical model. In ^ we extracted the weight factor uj from the experimental 
data. Now we have two input parameters for the analysis employing Models I-III. 
[see Eqs. (|2-l() - (|2-4() ]. In Models I and II, there is one model parameter (ej for Model 
I and an for Model II) corresponding to the two input quantities. A*™*' and w. We 
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determined the model parameters by minimizing 

rxtruc \ true] 2 r i2 

^ ~ (5At-)2 + (<5a;exp)2 ' ^ ' 

where A*™^ and Wexp were extracted from the experimental data in ^ and 0] The 
quantities (^Ag^p^ and Jwexp are the errors on A*™*' and Wexp; which are given in 
Tables IIIII and IIVI respectively. In Models I and II, X^^^^ and Ucai are functions 
of El and an, respectively, calculated using Eqs. (|2-1() and ()2-2() . In Model III, we 
solve Eqs. (|2-3|) and (|2-4j) for the given sets of Ag^p*^ and Wgxp to obtain eni and am. 
However, solutions of these equation may exist in unphysical parameter regions, such 
as em < 0. In such cases, we determine a "Best fit" solution by minimizing the above 
within the physical model parameter region, < em < 1 and < am < oo. 

5.1. Partially coherent model (Model I) 

We plot £i determined from Eq. H5-l|) as a function of the multiplicities of various 
collision systems in Fig. \^ We consider three cases: The open squares represent the 
results obtained from Eq. H5-l() using (A*''"'^, w), the solid triangles represent those 
obtained using {Xp^^,u;), and the closed circles represent those calculated from uj 
only, using the second equation in Eq. (|2-lj) . 

While the result denoted "a; only" exhibits clear increase of £i with multiplicity, 
which is consistent with the result obtained by the STAR,^^) neither (A*™°,cj) nor 
(Ap™'^,a;) exhibit such a clear dependence. This is because A*'^"'' in Table HTll does not 
display a clear multiplicity dependence and (JA*™^ <C 5wexp- The value of ei which 
minimizes is mostly determined by A*"^"^ and Ap™*^ only, i.e., the values of lo are 
not clearly reflected in the minimization of due to the fact that they possess larger 
errors than A*™*^. If this model is good enough, and if the experimental background 
for the correlation, such as the Coulomb correction is successfully removed, the result 
of (Ap™°,u;) should agree well with the "a; only" result. In Fig. |S1 these results seem 
to be quantitatively consistent. If the experimental accuracy of the three-particle 
correlation measures improves, these results will be more conclusive. 

5.2. Multicoherent model (Model II) 



Table V. The parameter qh in the multicoherent model (Model II) 



System 


From (A"'"°,a;) 


From (A»r,^) 


From ui only 


SPS S+Pb 




^■^'i-o.ai 




SPS Pb+Pb (NA44) 


40.121^3.84 


q rr+0.67 


' ^-0.38 


SPS Pb+Pb (WA98) 
RHIC Au+Au 


201.93+5^9 2 
10.541^1,0^ 


*.Oy_0.92 

2.86lJ;?| 


7 i:q+9-66 
'•Oy_6.53 



The results for an in the multicoherent model [Eq. H2-2|) ] are displayed in Ta- 
bled As in the case of Model I, three cases are shown. The results for (A*''"^, a;) and 
{Xp^°,uj) cases mainly reflect their A*''"^ and Ap"*' values due to the fact that 5A*''"° 
is much smaller than 6uj. In the (A^'''^'^, u;) case, an takes a very large value, coming 
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Fig. 5. £i as a function of the multiplicity. For the Pb+Pb collisions at the SPS, we plot both 
the NA44 and WA98 data with a slightly shifted horizontal axis for a clear comparison of the 
results. 



1 



from and to 
from A,pc^nd co 
from CO only 



from A 



true 



1. The value of an in the (A 



pc 



case is significantly smaher, but 



there is no clear multiplicity dependence in this case as expected from the fact that 
Ap",."*^ does not possess a clear multiplicity dependence. In the "w only" case, there 
are no solutions for the SPS S+Pb and Pb+Pb (WA98) data (the blank entries in 
Table IV|) . because this model has no corresponding value of an below lo ~ 0.82 [see 
Fig. 2 of Ref. 14)]. This implies that this model is not suitable for studying multi- 
plicity dependence; i.e., models with a chaotic background give a better description 
of the data. 

5.3. Partially multicoherent model (Model III) 

In the analysis using Model III, there are two output parameters (em and am) 
corresponding to the two inputs, A*^"*^ and w. In the following, we display the allowed 
regions of em and am which correspond to the sets of A*'^"^ and uj in Fig. El and Ap™^ 
and u! in Fig. [3 

In each figure, the lightly shaded area labeled "A" is the allowed parameter 
region corresponding to the value of A*''"'' or Ap"*', whose boundary is indicated by 
the solid line. Area B, bounded by the dashed (upper limit of uj) and dotted (lower 
limit of uj) curves, is the allowed parameter region corresponding to uj. The darkest 
area, which is the region in which Areas A and B overlap, is Area C, which is the 
allowed parameter region for am and em- The best fit values calculated from the 
values in Tables ITTTI and HVl are indicated by the squares. 

Figure ini^a) plots the result for the S+Pb collisions using A*'^"". It is seen that 
Area C is narrow, and the chaotic fraction has a lower bound near em = 0.6. The 
best fit value is em — 1 which reflects that A*'''^'' is near unity and a lower value of 
UJ. If we adopt a partially Coulomb corrected A*™'^, the situation changes. Due to 
the smaller Ap"^[Fig. [7|^a)], em is allowed for a wider region and it has a maximum. 
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The best fit value also shifts to Em = 0.67 and ani = 0.21. Note that, as we can see 
from Fig. [3 the upper bound of em is mainly dominated by the lower bound of to. 

Figure ini^b) displays the allowed region from A*"^"*^ and uj for the NA44 Pb+Pb 
collision dataset. Because the lower bound of A*''"° has a large value, 0.95, Area 
C allows both a mostly chaotic source with a small number of coherent sources 
and a number of the coherent sources with a small chaotic background. In this 
case, a solution of Eqs. (|2v-ij) - (|2-4j) exists in the unphysical region, em = 1.01 and 
am = 0.003. Hence, we determine the "Best fit" point by minimizing in Eq- (|5-1() . 
The result of the minimization gives the "Best fit" at em = and am = 40.1. This 
result does not always imply that the multicoherent picture is good; the difference 
between the value of for this minimum and in another case, for example, em = 0.99 
and a = 0.01, is much smaller than unity. 

If we adopt the partially Coulomb corrected chaoticity Xp^^, Area C in Fig. [7|^b) 
becomes narrow, as in the S+Pb case. Because the lower bound of uj is larger than 
that in the S+Pb case, the upper limit of em becomes smaller. Similarly, this case 
does not have a solution of em and am within the physical region (em = 1.56, am = 
1.12), and therefore the "Best fit" point is determined by minimizing in Eq. (|5-1|) . 
Though the location of the "Best fit" point corresponds to the multicoherent picture 




D 0.2 OA 0,6 tt.'A 0.2 OA 0,6 O.-A 

Fig. 6. Allowed regions for A*"^"" in the various systems. (See text for details.) 
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Fig. 7. Allowed region for Xpc° in the various systems. (See text for details.) 



(em = 0), the opposite case (em ~ 1) is statistically allowed for the same reason as 
in the previous case. 

Similarly, Figs.lHl^c) andEt^c) display the results for the WA98 dataset. Since the 
central value of A*'''^'^ exceeds unity, only a very small region near eni = 1 is allowed 
in Fig. Efc). This situation changes drastically if one adopts the Ap™*^. In Fig. Etc), 
the allowed region has a shape similar to that in the NA44 case [Fig. Etb)] . The best 
fit value is located at eui = 0.47 and am = 1-29. 

Finally, we display the results for Au+Au collisions at RHIC in Figs. Efd) and 
IT^d). From A'^'^"'^ and lo, it is difficult to distinguish the structure of the source; 
both a large chaotic fraction with a small number of coherent sources {em ~ 1 and 
o^iu < 1) and small chaotic fraction with a large number of coherent sources (em ~ 
and am > 1) can reproduce the experimental data. The best fit value is located 
at em = 0.75 and am = 0.77. However, if we use the partially Coulomb corrected 
two-particle correlation data, the allowed region is strongly restricted. Although 
the solution of Eqs. H2-3() - H2-4|) exists in the unphysical region, em = —0.54 and 
am = 6.83, the "Best fit" point is located at em = and am = 2.86. The maximum 
value allowed for em is 0.4. This means that the strong chaotic behavior observed 
at the RHIC is due to the production of a cluster of coherent sources. 
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As shown in the Fig. the maximum bounds on eni are mainly determined 
by the lower bounds on w, which become larger as the energy and the multiplicity 
increase. The fact that am becomes larger as the multiplicity increases is consistent 
with the result of a previous analysis by one of the authors (H. N.), in which the 
chaoticity and the weight factor were experimentally determined. ^^-^ In Fig. |H1 we 
plot the maximum and minimum values of am for the sets of X^^^ and u; extracted 
from Fig. d We can see that minimum value of am increases with the multiplicity. 
This suggests that, as the collision energy and the multiplicity increase, the number 
of coherent sources increases. Note that the maximum values are the same as those 
obtained using Model II for X^^^ (see Table IVj) . because am takes its maximum 
value when em = 0. Furthermore, this tendency seems to be correlated with the 
plausibility of the statistical model (x^/A^dof)) given in Tables Q] and ^ Though 
we do not know any explicit relation between Models I-III and the hadronization 
mechanism, this result may reflect a possible hadronization mechanism from the 
quark-gluon plasma phase created in the collisions. 

§6. Summary 

In summary, we have investigated the degree to which the pion sources are 
chaotic in various heavy ion collisions by analyzing the two- and three-particle cor- 
relation data with three kinds of particle production models. For the two-particle 
correlation data, we have extracted the "true" chaoticity by considering long-lived 
resonance contributions to the pion multiplicity, with the help of a statistical model. 
Using simple source functions, we simultaneously investigated the two- and three- 
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Fig. 8. The minimum (solid squares) and maximum (solid circles) value of am in Model III for 
various collisions as a function of the multiplicity. For the Pb+Pb collisions at the SPS, we plot 
both the NA44 and WA98 data with a slightly shifted horizontal axis for a clear comparison of 
the results. 
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particle correlation functions to extract the weight factor u; of the three-particle 
correlator. Incorporating the chaoticity and the weight factor into the models, we 
have studied the chaotic fraction and mean number of the coherent sources. The 
results for ej obtained from uj indicates that the system becomes chaotic as the mul- 
tiplicity increases. This result is consistent with Ref. 51). From a multicoherent 
source point of view, it is concluded that pions at higher collision energies may be 
emitted from a cluster of coherent sources and the number of sources increases as 
the collision energy and the multiplicity increase (Fig. |HJ . 
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